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^ ■ Abstract 



The Gas Pixel Detector was designed and built as a focal plane instrument 
for X-ray polarimetry of celestial sources, the last unexplored subtopics of 
X-ray astronomy. It promises to perform detailed and sensitive measure- 
ments resolving extended sources and detecting polarization in faint sources 
in crowded fields at the focus of telescopes of good angular resolution. Its 
polarimetric and spectral capability were already studied in earlier works. 
Here we investigate for the first time, with both laboratory measurements 
and Monte Carlo simulations, its imaging properties to confirm its unique ca- 
pability to carry out imaging spectral-polarimetry in future X-ray missions. 
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1. Introduction 

The Gas Pixel Detector (GPD) [H, 0, i has been designed and built to 
perform time resolved imaging spectral-polarimetry of X-ray celestial sources 
by means of the photoelectric effect. It is the longsought quantum leap in 
sensitivity with respect to the classical techniques [4j . Successively the same 
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method was implemented exploiting an alternative technique without any 
2-D imaging capability Jj] • 



While the polarimetric and the spectroscopic performances of the GPD 
have been already studied in detail ja, [^j , its imaging properties have been 
estimated, so far, only by Monte Carlo simulations evaluating the point of 
conversion of the impinging photon with a suitable algorithm j^]. Still the 
GPD was devised primarily to be exploited at the focus of conventional and 



multi-layer hard X-ray optics [22j. At this pur pos e the GPD has been consid 



ered for different proposed space missions [14], ll2l, |13|, ll5j . While X-ray optics 



with the exquisite quality of those of Chandra will not be again available in 
the near future, for POLARIX (or XIPE, proposed in 2012 to ESA as a small 
mission for a launch on 2017) the JET-X optics (three complete mirrors, two 
Flight Models, FM, the heritage of the former Spectrum X-Gamma mission, 
and one Qualification Model), can still be used. They have a Point Spread 
Function (PSF) with a measured Half Energy Width (HEW) at 1.5-keV of 
15.1" (FM-1) and 14.6" (FM-2) Q. 

We report on the first measurement, with a narrow beam at 4.5 keV, of 
the position resolution of a GPD filled with a mixture of He-DME 80-20 (1 
bar) with a large Gas Electron Multiplier (GEM), This proves the capability 
of the GPD to perform space-resolved measurements, such as those of Pulsar- 
Wind Nebulae, Supernova Remnants or X-ray jets with a suitable accuracy, 
and to observe the faintest source accessible to X-ray polarimetry especially 
in crowded fields 2(|. We also compared such measurement with the Monte 



Carlo simulation relating the characteristics of the image to the modulation 
factor. Finally we showed the effect on imaging of the non-uniformity of the 
drift electric field in case the GEM plane, and in particularly the guard-ring, 
is not powered. 



2. The GPD prototype and the relevant algorithms 

The GPD is a gas cell with a beryllium thin window, a drift region, a 
perforated multiplication plane (the GEM) and a collecting plane with an 
hexagonal pattern. 

The collecting (or readout) plane is the top layer of an ASIC CMOS 
with a square surface of 1.5 cm x 1.5 cm having 105600 hexagonal pads of 
50 /im pitch. It is designed with a very high precision due to the micro- 
electronics technologies employed. The positions of the hexagonal pads set 
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up, therefore, a very high accurate reference system which is essential for an 
accurate imaging. 

Hereafter we summarize, to understand the considerations in the following 
sections, the sequential steps of the algorithm that processes the signal from 
each hexagonal pad that collects the charges of the track. The algorithm is 
aimed to evaluate: 



• the barycentre using all the pixels; 

• the angle, for which the second moment, with respect to the barycentre, 
is maximum; 

• the third moment (skewness), also with respect to the barycentre, to 
determine which of the two end points contains the impact point (by 
the smaller charge density); 

• the impact point, using the pixel only within a sub-region of the initial 
part of the track as it is identified by the skewness. The extension of 
this sub-region depends on the gas mixture and it is evaluated by using 
the maximum second moment as a weight; 

• the angle (emission direction) that maximizes the second moment with 
respect to the impact point weighting the charges according to their 
distance from it. 

The first sealed GPD [i| had a GEM plane with the perforated surface 
as large as the ASIC CMOS chip and a guard-ring, not perforated, just 
slightly larger. The present prototype 1^| has, instead, a much larger guard- 



ring extending as far as 2.25 cm from the edge of the ASIC chip. This latter 



design, see fig. 1(a) and fig. 1(b) for a comparison, allows for a better control 
of the uniformity of the electric field. Border effects, other than those due 
to the extension of the track and to its incomplete imaging close to the edge 
of the sensitive region, are now therefore expected to be vanishing. We will 
show later how the guard-ring influences the position resolution of the GPD. 

The GPD has a drift/absorption length of 1 cm and a beryllium window 
50 /iin thick. The GEM employed with this new design is 50 /xm thick, with 
holes of 30 /iin of diameter with 50 /im pitch. The holes are laser etched on 
a liquid crystal polymer substrate. The GEM is manufactured by SciEnergy 



under the supervision of RIKEN [16 
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(a) 



(b) 



Figure 1: (a). The old design of the GPD with a narrow guard-ring in the GEM plane 
(b). The present design with a larger guard-ring (l9| . 

3. The X-ray beam profile 

In order to measure the position resolution we arranged a narrow X- 
ray beam by means of a collimator system connected to the exit beryllium 
window of a 50 W (50 kV, 1 mA) commercial X-ray tube with a Titanium 
anode by Oxford Instruments (series 5000) having an anode spot smaller 
than 150 /im. The X-rays from the tube, operated at 30 kV and 400 //A, 
were collimated by means of two diaphragms, one made of brass having a 
diameter of 500 /im connected to the exit window of the X-ray tube and 
one made of tungsten with a diameter of 25 /im on the opposite side close 
to the detector window. The two diaphragms were interconnected by two 
coaxial brass shielding tubes and separated by 17 cm. The large distance 
necessary to reduce the beam divergency does not allow for measuring the 
position resolution at energies below the Titanium X-ray K lines at 4.51 keV 
(K a ) due to the intervening air absorption. In order to measure the actual 
beam size, we scanned the X-ray beam with a brass slit with sharp edges, 
at the same distance of the GPD, in two orthogonal directions and at a 
constant speed using as detector a commercial one-pixel Si-PIN XR-100-CR 
manufactured by Amptek. The geometrical surface of the detector is 5 mm 
x 5 mm. This system, when interfaced to the two orthogonal linear stages 
with 0.5 /im resolution (model ILS50CC), manufactured by Newport and 
operated via computer by means of an 8-axis XPS controller |2l|, allows for 
continuously shifting the detector and the slit with respect to the beam at a 
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minimum speed of 1 /im s . Data from the Si-PIN detector were acquired 
in rate-meter mode with a time-bin of 1 s, selecting a window in energy to 
include only the titanium K-lines. With this set-up, the time needed for 
obscuring the beam at the slit passage provides the beam size along the 
direction of the scan. 

If the beam is described by a 2-D gaussian, as in the following equation 

H): 

, (X-X ) 2 (y-y ) 2 - 

f(X,Y) = Ae (1) 

the counting rate as a function of time (or as a function of position being 
the speed constant and known) resulting from the scanning along the X (or 
Y) direction can be described by the erf function as in eq. EJ This is, 
actually, the primitive of eq. [H after the integration between -oo and +oo of 
the counts in Y (or X) direction. 

The counting rate during the scanning in Y and X are shown in figure [2] 
together with the fitting function : 

!(X) = ^[1 - ern ^-X'k (2) 

Such fitting procedure allows for deriving the relevant parameters (ax and 
cry) of the X-ray beam and comparing them with the results of the section 
HJ The results of the fit are shown in the table [TJ 

Table 1: Result of fitting procedure of the beam scan. The xt 1S 1-0 for the X scan and 
1.12 for Y scan. 



Param 


Fit results 


A x 


(171.71 ± 1.21) Counts/s 


X 


(132.42 ± 0.31) fim 


(?x 


(8.7 ± 0.3) /im 


A Y 


(171.67 ± 1.16) Counts/s 


Y 


(152.41 ± 0.78) /im 


a Y 


(14.7 ± 0.4) fjm 



The beam shape is therefore consistent with being gaussian with two dif- 
ferent elongation in the two orthogonal directions. Such different elongations 
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may be due either to the position of the 500 fim diaphragm with respect 
to the anode spot of the X-ray tube or to a possible residual misalignment 
between the two small diaphragms at the opposite end of the two collimator 
brass tubes. 




100 200 300 100 200 300 

Seconds (1/2m/s) Seconds (I/im/s) 

(a) (b) 

Figure 2: (a). Scan along the Y axis. The result from fitting eq[5]is ay = (14.7 ± 0.4) 
/im (b). Scan along the X axis. The result from fitting eq[5]is ax = (8.7 ± 0.3) /zm. 

We did not investigate further the cause of the two different sizes since 
their values are already smaller than the expected position resolution of the 
GPD. 

4. The position resolution determination 

To study the imaging capability of the GPD we made a set of three 
measurements shifting it by 300 /im with respect to the central position 
(the center of the ASIC chip readout plane) in orthogonal directions. This 
central position was, actually, found manually by iterating once the shift and 
the derived centroid of the image of the beam. 

The counting rate during each acquisitions was about 10 counts s _1 , 
by far larger than the background counting rate, for a total of about 8700 
counts/acq., all windowed in energy to include only the Ti K photo-peak. 
The images of the three acquisitions all together are displayed in figure [3] 
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Figure 3: The real images of the beam detected by the GPD. These images were ob- 
tained by shifting the GPD by 300 /jm in orthogonal directions with respect to the central 
position. 



where they appear well separated. A closer view shows that they are char- 
acterized by a dominating core with external wings. We investigated the 
properties of the core and the wings in detail later in section [5] by means of 
Monte Carlo simulations. In this section we, instead, measured the proper- 
ties of the core. A bi- dimensional gaussian including two free ax and ax 
parameters (see eqJI]) was fitted to the data considering only the central part 
of the impact point distribution in a frame of 120 /im x 170 /iin around the 
peak position. The fit was performed by using the MINUIT package included 
in the ROOT 5.0 environment. The results of the fit are shown in the table 
|2j The x 2 °f th e fit of the core was 200.9 with 196 degree of freedom. A 
close-up view of the image of beam in the central position and the fitted 



function are shown in figure 4(a) and 4(b) 



Having determined in section [3] the main parameters of the beam pro- 
file we are in the position to correct the measured values. We subtracted 
quadratically the size of the beam as derived in table [T]from the position res- 
olution of the GPD as in table [2J The results for the two standard deviations 
are a x PD = (28.1 ± 0.5) and ay FD = (34.3 ± 0.6) fim. The remaining 
significative, albeit small, difference in size between the two directions could 
be due to the different sampling in X and Y because of the hexagonal readout 
pattern. While X is sampled with 50 /im resolution, thus the hexagonal pitch, 



Table 2: Fit results of the real data acquired by the GPD in the central position, in a 
frame of 120 /im x 170 /jm around the peak position, using the 2-D gaussian of equation 

m 



Param. 


Fit results 


A 


95.36 ± 1.91 counts 


x 


(2.1 ± 0.4) /im 


Yo 


(-17.4 ± 0.5) /im 


ox 


(29.5 ± 0.4) /im 


a Y 


(37.3 ± 0.4) /im 



Y is sampled with a resolution of 50 /im x cos(30°) = 43.3 /im. However we 
cannot exclude, also, a residual misalignment of the X-ray beam axis with 
respect to the GPD plane and the effect on the different elongation due to 
the inclined penetration. The control of such angle is, indeed, mostly relied 
on machined parts and an alignment procedure based on measurements at 
controlled different inclinations will be studied. For the time-being we retain 
that <Jx PD an d Cy PD are close enough to assess that no major problems in 
the set-up, in the CMOS ASIC read-out or in the algorithm are present. 
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Figure 4: (a). A close-up of the image of the central beam. (b). 2-D Gaussian function 
fitted to the core of the image. 
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5. The Monte Carlo simulation of the beam profile 



A simulation was performed with the Monte Carlo we routinely use to 
anticipate the polarimetric performances of the GPD. The Monte Carlo sim- 
ulates the generation, the transport, the multiplication and the collection 
of the Auger-electrons tracks and of the photoelectron tracks from different 
atom shells. Such tracks are generated by X-rays of a given polarization 
degree and of a given energy in a given gas mixtures. The main routines of 
the Monte Carlo were earlier described in jijl, 0, [h| . The simulated data are 
processed by algorithms equivalent to those used for the real data collected 
by the GPD. The agreement between the Monte Carlo and the experimental 
data in terms of modulation factor was already studied for different gas mix- 
tures and detectors, see @, 0, 19]. Here, instead, we compare the position 
resolution. At this regards we simulated a narrow pencil beam of 4.5 keV un- 
polarized photons impinging on the center of the GPD. For each simulated 
track we derived the reconstructed impact point, together with the recon- 
structed emission direction. The total number of counts in the simulation 
is larger with respect to the acquired real data of a factor of about 3. The 
location of the impact points resulting from the simulation is shown in figure 
[5j The presence of a core and external wings, evident in this image, will be 
characterized in section 15.11 and 15.21 In this paragraph we show the results 
of the fitting procedure to determine the parameters of the simulated image 
profile. 

Such image profile, together with the 2-D gaussian function used to fit 
the data (within a square wide 110 fim x 110 /urn), is, indeed, displayed in 



figure 6(a) and 6(b) respectively. 



Table 3: Results of the fitting of the core of the simulated beam image with the 2-D 
gaussian function of equation [T] 



Param. 


Fit results 


A 


363.9 ± 4.4 counts 




(0.00 ± 0.02) fxm 


Vo 


(0.00 ± 0.08) Aim 


Ox 


(30.8 ± 0.3) ^m 


Oy 


(31.9 ± 0.3) ^m 



The results of the fit are in table [3] and the xl is now 2.2. Such a large 
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Figure 5: Scatter plot of the reconstructed impact points from simulated data. Wings are 
defined as the set of impact points outside a core of 200 /jm of diameter. 
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Figure 6: (a). Monte Carlo simulation of a narrow beam of 4.5 keV as imaged by the GPD. 
(b). Gaussian fit of the core (within a square wide 110 /im x 110 /Km) of the simulated 
data. 
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value may be is due to the contribution of the wings which are included in 
the corner of the selected region. Actually fitting the same function on a 
larger set of the data (a square wide of 130 /zm x 130 jum), the xl increases 
to 2.5 probably because of the larger contribution of the wings. 

If we compare the results of table [3] with cr^ PD and crp PjD , we see that 
they agree within about 3 /im. The agreement is very good also considering 
that the simulated beam is aligned by construction with respect to the GPD 
plane. 

5.1. The core and wings in the simulated image 

Puzzled by the presence of the wings, we used the Monte Carlo data to 
investigate their origin and to study their contribution to the sensitivity for 
polarimetry. With the Monte Carlo is, indeed, easier to relate the evolution of 
the original simulated track to the event detected by the GPD after transport, 
blurring, multiplication and collection. The exact definition of the core and 
the wings is beyond the scope of this paper. Therefore we define, somewhat 
arbitrarily, the core as the region inside 200 fim in diameter from the peak 
position, and the wings as the set of the impact points outside the core. 
With this definition in mind we compared the detected number of events 
with those expected by the 2-D symmetrical gaussian function. We recall 
that the density of probability (equation [T]) with the normalization constant 
27r ^ ff2 , if circularly symmetric, can be written as : 

1 _ r 2 

f( r ) = — x r x e ^ (3) 
a z 

The fraction of events outside a circle of radius ro can be written as in 
equation @] : 

roo 1 
r 



f(r)dr = e 2^2" (4) 

The fraction of the reconstructed impact points outside a radius of 100 /im 
in the Monte Carlo data is indeed 16.6 %, while the expected fraction is only 
0.5 %. This excess represents the contribution of the wings. The investigation 
of the wings is important because some wings are, generally, present also in 



the Lorentian PSF of the X-ray optics [17], |T8| and an interplay between the 



two could be expected in focal plane experiments. 
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An inspection of the Monte Carlo tracks with calculated positions outside 
the core showed that the presence of wings is due, indeed, to an incorrect de- 
termination of the reconstructed impact point. In practice the reconstructed 
impact point is on the opposite side with respect to real impact point. This 
happens because the algorithm calculates the skewness for determining the 
end point of the track (see section [2]) containing the Bragg peak that is the 
end of the photoelectron evolution. The distribution of charges for some kind 
of trajectories provides blurred collected tracks with the skewness giving the 
'wrong' sign. Consequently, instead to identify in the track the correct end 
point, the incorrect one is derived. A simulated track where this inversion 
happens is shown in figure 7(a) where it is displayed its projection on the X 
& Y plane (the GPD readout plane) in the Monte Carlo coordinate system. 
In figure 7(b) it is shown the image of the same track blurred by diffusion, 
multiplied by the GEM and collected by the readout plane. The green cross, 
in the latter, is the original impact point, the blue cross is the barycentre 
and, finally, the red cross is the reconstructed position. An inspection of 



the original track shows, see figure 7(a) , that the impact point in the original 



frame coordinates is at (X = 0; Y = 0), the Auger-electron evolves rightwards 
while the photoelectron firstly evolves in the same semi-plane containing the 
Auger-electron, secondly scatters on a nuclei and moves on the opposite semi- 
plane. The derived distribution of the blurred charge is such that a larger 
density is at the end point containing the original impact point and not at 
the opposite end point containing the Bragg peak: the evaluation of the 
skewness, therefore, produces an inversion. A consequence of such an error 
is that the emission direction is, also, changed (see again figure 7(b) and its 
caption) and almost inverted. Due to the 180° symmetry a certain, albeit 
smaller, sensitivity to polarization is expected to be preserved also in the 
events contained in the wings. This is the object of the study presented in 
the section 15721 



5.2. The core and wings: the impact on the polarization sensitivity 

We studied also by Monte Carlo simulations the contribution of the core 
and wings to the polarimetric sensitivity of the GPD. Such study, in our lab- 
oratory, at the present time can only be done by simulation. The generation 
of polarized radiation, in our facility, greatly reduces the flux from the X-ray 
tube and, therefore, the use of a collimator system with a narrow beam in 
addition to the polarizer is of no practical use. 
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Figure 7: (a). A simulated track in the wings projected on the GPD plane, (b). The 
same simulated track after transport, blurring, multiplication and collection. The green 
cross and the arrow line is the original impact point with the emission direction while the 
red cross and arrow line is the reconstructed impact point and the reconstructed emission 
direction. 



However the agreement between the Monte Carlo data and the experi- 
mental data, both in terms of modulation factor and position resolution, is 
such that a comparative study of the wings and the core in terms of polar- 
ization sensitivity is reasonable in this way. 

We simulated for this study a zero-width narrow X-ray beam of 4.5 keV 
and 100 % polarized and we derived the reconstructed impact point and the 
emission direction for each event. The histogram of the emission directions, 



the so called modulation curve, is reported in figure 8(a) and in figure 8(b) 



for the core and the wings respectively with the following function to derive 
the modulation factor : 



/(0)=A + iW(0-</> o ) (5) 

where A and B and 4>o are parameters of the fit. The modulation factor 
fi is given by the equation [6] : 



B 



// 



2A + B 

The modulation factor of the whole data set is /i 
contribution of the core is /i = (55.3 ± 0.6) 



(6) 

(50.9 ± 0.6) %. The 
o and the contribution of the 
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Figure 8: (a). Modulation curve from simulated 100% polarized 4.5-keV photons in the 
core. The modulation factor is fJ, = (55.3 ± 0.6) %. (b). Modulation curve from simulated 
100% polarized 4.5-keV photons in the wings. The modulation factor is [i — (29.6 ± f .6) 
%. 



wings decreases down to fi = (29.6 ± 1.6) %. We therefore demonstrated 
that the core, as expected, has a larger polarization sensitivity with respect 
to the wings. This is interesting because, by selecting the events contained in 
the core, we could determine the polarization of a point-like celestial source 
with a better confidence. We recall, actually, that the polarization sensitivity 
is proportional to the so called quality factor (equation [7j) : 

QF = fjLy/i (7) 

where e is the efficiency. In our case QF core = 0.505 and QF a u = 0.509 
therefore basically identical, however the larger modulation factor obtainable 
selecting only the events in the core allows for a slightly better control of the 
incidence on the measurement of possible residual systematics. Moreover 
the lower, but not zero, modulation factor of the wings allow for using these 
events once properly weighted. In any case the analysis shows that the main 
origin of the wings is that, for a fraction of events, mainly those characterized 
by photoelectrons that underwent to a large scattering angle with an atomic 
nucleus, the skewness is not sufficient to identify the section of the track that 
includes the impact point. This means that an improved pattern recognition 
can be mitigates significantly the problem. 

The effects on sensitivity of the selection of the core and the wings needs a 
more accurate study. For example they are expected to be energy dependent 
through the different average size and shape of the tracks. Such study is 
beyond the scope of this work, but we recall that spectral capability of the 
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Figure 9: Distortion of the two images of the beam shifted 300 fim after the disconnection 
of the upper guard-ring. 

GPD allows for such an investigation. 

6. The effect of the large guard-ring on the image profile 

The two different configurations up to now used for the GPD, the old 
configuration, with a small guard-ring in the GEM plane, and the present 
configuration, with a larger guard-ring, were already described in section 
EJ We can not directly compare the two configurations in terms of position 
resolution, since the old configuration is not available at the moment, but 
we can show the effect of the presence of the guard-ring on the location 
capability of the detector. Therefore, to study how the uniformity of the 
electric field impacts on the image of the narrow beam, we disconnected the 
upper guard-ring of the GEM plane. The results are shown in fig. [9] for two 
of the three positions of fig. [31 The effect is indeed enormous. The image of 
the beam appear now larger, distorted and rotated even at the center of the 
active region at about 7.5 mm from its border. 

The electric-field, with this disconnection, is not the same of that of 
the old-one, with a smaller guard-ring and a smaller aspect ratio (the ratio 
between the thickness of the drift region and the linear dimension of the 
surface area of the GEM plane). However this measurement suggests that 
the present configuration assures better performances in terms of drift-field 
uniformity and position resolution capability. 
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7. Conclusions 

We studied the position resolution of the GPD by means of laboratory 
measurements and Monte Carlo simulations, at 4.5 keV, using a narrow X- 
ray beam. After subtracting the contribution of the beam size which was 
measured by means of a scanning technique, we obtained a position resolution 
more than twice better than that already found with a Micro Pattern Gas 
Detector with 200 /im pitch and a GEM with a larger pitch [9j . The measured 
position resolution is very close to the Monte Carlo results. By using the 
latter and equation HJ the Half Energy Width (HEW, spatial) found is 36.3 
/im. The expected HEW (spatial) by the JET-X optics (with a focal length of 
3500 mm and a HEW (angular) of 19.3" [8]), is 327 /im, a factor of nine larger, 
including also the effects of inclined penetration of the photons collected in 
the absorption depth of the GPD. The angular resolution of an experiment 
with the GPD at the focus of such X-ray optics is, therefore, driven primarily 
by the quality of the latter. 

We detected the presence of wings in the image of the beam in both the 
real data and the Monte Carlo data. Such wings were studied primarily by 
simulations and are due to an inversion in the location of the impact point 
due to peculiar event trajectories and to the blurring. The core events have 
a larger modulation factor with respect to the events in the wings suggesting 
that a spatial selection of the data, a different weight, or a better pattern 
recognition, especially for polarimetry of point-like celestial sources, would 
lead to an improved statistical significance. 

Finally we showed the effect of the presence of the large guard-ring on 
the image profile of the narrow beam. 
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